z1 <- (61.4 - 60)/4
z2 <- (63.5 - 60)/4

prob <- abs(pnorm(z2) - pnorm(z1))
prob


#x <- rnorm(10000, mean=60,sd=4)
#hist(x)